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FOREWORD 


The research reported herein was performed in the Materials 
Science Department of the University of Virginia, Charlottesville, 
Virginia 22901. The work was initiated under grant number 
NASA-NAG-1-419, "A Study of the Applicability of Nucleation Theory to 
Quasi-Thermodynamic Transitions of Second and Higher Ehrenfest- 
Order." This report covers research performed from December 1983 to 
December 1984. Since the major portion of the first year's work 
involved the theoretical development of a relatively unexplored area the 
results of these efforts are presented in a more detailed fashion than 
usual in this Annual Report. The investigators for this project are 
Dr. R. Edward Barker, Jr., principal investigator, and Kenny W. 
Campbell, graduate research assistant. 
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SECTION I 


ABSTRACT 

The applicability of classical nucleation theory to second (and 
higher) order thermodynamic transitions in the Ehrenfest sense has 
been investigated and expressions have been derived upon which the 
qualitative and quantitative success of the basic approach must 
ultimately depend. The expressions describe the effect of temperature 
undercooling, hydrostatic pressure, and tensile stress upon the critical 
parameters, the critical nucleus size, and critical free energy barrier, 
for nucleation in a thermodynamic transition of any general order. 
These expressions are then specialized for the case of first and second 
order transitions. The expressions for the case of undercooling are 
then used in conjunction with literature data to estimate values for the 
critical quantities in a system undergoing a pseudo-second order 
transition (the glass transition in polystyrene). Methods of estimating 
the interfacial energy Z in systems undergoing a first and second order 
transition are also discussed. Plans for future research and a list of 
publications and presentations accompanied by manuscripts and 
abstracts are also included. 


1 


SECTION II 


INTRODUCTION 


A. General Background 

As in the case of so many other important concepts in science, the 
foundations of nucleation theory were laid a century ago by Gibbs. ^ 
Important developments which put the subject of homogeneous nucleation 
in a form recognizable by modern students occurred in the 1920's and 
1930's when Volmer and Weber^ and Becker and Doring^ sought to 
understand the phenomenon of vapor condensation. Subsequently the 
concepts of heterogeneous nucleation were developed, so that modern 
nucleation theory includes both homogeneous and heterogeneous models. 

The theory has been applied to phase transformations in metals 
and to a large extent metallurgists have tended to dominate in the major 
developments in the field for many years. Nevertheless, some signifi- 
cant specialized applications of nucleation theory have been made by 
polymer scientists. Among the most notable are the efforts by 
Price^, Hoffman and Lauritzen^, Lindenmeyer^, and others to 
understand the originally unexpected lamellar morphology of polymer 
single crystals and the kinetics of chain folding in such crystals. 

Nucleation theory has previously been successfully applied only to 
thermodynamic transitions of first order in the Ehrenfest sense. 
However, in this work, classical nucleation theory is applied to 
transitions of second (and higher) order. 

B. Meaning and Significance of a Second Order Transition 

In the Ehrenfest classification of phase transitions, the order of a 
given transition is determined by which derivatives of the Gibbs free 
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energy are discontinuous. Thus, in the present thermodynamic con- 
text, the phrase "seccnd order" refers to physical processes of a 
fundamental nature and not necessarily to correction factors of small 
magnitude. According to the definition, a second order phrase transi- 
tion in the Ehrenfest sense is one in which the Gibbs free energy and 
its first derivatives are continuous at the transition point while the 
second derivatives are discontinuous. The primary discontinuous 
properties associated with the second derivatives are changes in the 
thermal expansivity Aa, the isothermal compressibility AfJ, and the heat 
capacity at constant pressure AC . Hence, a second order phase 
transition is characterized by "jumps" in such thermodynamic quantities 
as those mentioned above and thus differ from a first order transition 
only in the quantities which exhibit discontinuity and not in fundamen- 
tal significance as thermodynamic phase transition. 

C. The Thermodynamic Basis of Nucleation Theory 

As is widely known, when a pure system of molecules (atoms 
included) which exists as an equilibrium (stable) phase a at a given set 
of intensive variables (T^, p^, etc.), is subjected to new intensive 
conditions (T 2 , pj, etc.) where a is no longer the only stable phase 
but where one or more new phases are stable, then it is possible that 
the expected transformation 

a- 3 (1) 

will occur only slowly, if at all, within a finite observation period. 
According to the precepts of nucleation theory, the reason for this 
inhibition of the expected transformation is that molecular fluctuations 
and diffusive motion must occur which create a small particle of 3 within 
the e-milieu (but see ref. 7) and :hat there is an interface between the 
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a and 8 regions in which the molecules have an excess free energy V * 
per unit area. Thus the true Gibbs free energy change 
corresponding to the transformation of a small amount of a into a small 
nucleus of 0 is not merely the so-called bulk thermodynamic value AG 
(energy per mole); it is instead 

AG - v e (AG/u p ) + A a(J ar ap (2) 

where is the volume of the 0-nucleus, is the molar volume, and 
is the interfacial area of 0. If the growing nucleus has more than 
one type of interface, then the term A c would imply a summation 

op op 

over the relevant areas. 

Thus, according to nucleation theory local density fluctuations will 
lead to the formation of a small nucleus of 0-phase within an a-phase 

natrix if the temperature is less than the normal transition temperature 

o o 

T^ for a ^ ivhere 6 is the stable phase when T < T^ (Fig. 1). The 

total intertacial energy A^JT^ will oppose the enlargement of the 

nucleus, and a bulk free energy gV^ will encourage growth (Fig. 2). 

The term 

* 1 <Gj - G„)/o 6 (3) 

represents the free energy change per unit volume of transformed 
material. Equation (2) now becomes 

AG * gV p + TAg . (4) 

D. The Customary Approximations 

In the standard treatments, AG s AH - TAS is approximated, at a 
degree of under cooling 

0 * T° - T (5) 

o o o 

by AH * AH^, AS * AH^/T^ which may be seen to give 
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• « 

AG * AS t • 0 where AS t 

or 

0 

g * * t 0 with s t = AS t /Up. (7) 

The standard treatments next proceded by taking 3(AG)/3r = 0, where 
r is the radius of a spherical nucleus for which 

Vjj * (4/3)xr 3 and = 4sr* (8) 

This gives the size of the critical nucleus r fi to be 

r c » “21/g . (9) 

For the stable a ■* 0 transition, g will be negative (g s 0 for the 
equilibrium transition). Substitution of Eq. (9) back into Eq. (4) gives 
the usual type of expression for the "activation barrier" that must be 
overcome before the continued growth of the 0-nucleus can take place. 

AG C * (16Tr/3)I 3 /g* * 16.76 r 3 /g 2 . (10) 

This barrier, and the critical size r , are of course both dependent on 

c 

the degree of under cooling, through Eq. (7), therefore 

-2 r(T/ 0) 

r c = -2T/s t 0 » jji and AG C - <*§*) *'/ s*0* . (11) 
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SECTION III 


SUMMARY OF PROGRESS, PAPERS AND PRESENTATIONS 

A. Summary of Progress on Research Grant NASA-NAG-1-419 

Tha main idaa behind the basis for our work is that if we break 
away from the conventional approximation that leads to Eqs. (6) and 
(7), instead introducing the assumption that AG and therefore g can be 
expanded as a Taylor series in 6, to which Maxwell's relations and some 
other thermodynamic operations are subsequently applied, then we have 
a powerful, relatively simple, and (as far as we have been able to 
ascertain) a previously undiscovered formalism for predicting (1) 
non-linear correction terms for first order thermodynamic transitions, 
(2) the size and barrier height for second order transitions (in the 
Ehrenfest sense), and (3) the possibility of a formal extension of 
nucleation theory to even higher order transitions. Since the new 
theory brings in the higher order effects as analytic thermodynamic 
terms (rather than merely mathematical or empirical approximations) it is 
a more flexible theory for predicting other effects, such as the influ- 
ence of changes in pressure, stress, electric field, etc. The new 
theory, for the same reason, should also be more amenable to the 
incorporation of ideas from statistical mechanics; and to compositional 
and kinetic extensions such as Couchman's 1 * T theory and the 

Sr 

(9) 

time-temperature transformation concepts of Enns and Gillham 
1. Expansion in Terms of Undercooling. 

The initial situation to which the formalism described above is 
applied is a system which is temperature undercooled below a normal 
(STP) transition point. 
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(a) Centra! features of the new approach. 

As outlined in the previous section it is possible to 
construct a useful, and apparently overlooked, formalism by what in its 
barest 'orm is the expedient of expanding the effective free energy of 
the transformed nucleus in a Maclaurin series in 0 (the degree of 
undercooling), and then utilizing thermodynamical relations to evaluate 
the expansion coefficients. For example, if we take 

8 * 8(0) * 8(0) + 8’ (0)0 + (l/2)g"(O)0* + . . . (12) 

* a Q + a x e + a 2 $* ... (13) 

then g’(0), which means (ag/30) p at 0 - 0, can be found as follows. 


but 


so that 


where 


'■ae ; * 

p p p p 


e 6 ’ G . * H S • H . ‘ T ‘ s > - «.)• 


-* 1 - <v' v e " h ' Te ’ 


h 1 (H p - H a )/v, and a 5 (S, - SJ/v,. 
Carrying through the operations required. 


8 — - j ( . ( “) 

v fl 1 3T J '‘3T ; , 

P & L P Pj 


1 

t 

»G i 

a. 


1 3v 6 

- V? <V V <«*> 

p p 

“ ^ ( V s « ) ■ «v 


(14) 

(15) 

(16) 
(17) 


(18) 

(19) 
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where ut* hat been mada of tha tharmodynamic ralationt S * -OG/iT) p 
and tharmal expansivity * a * (*v/3T) p /v. 


(b) Calculations of tha coefficients. 
Further reduction of Eq. (18) gives 


(jj) * * t * ( 5t^ 

P P 

and thus, using Eq. (13) at T * T^, or 0 * 0, we obtain 


( 20 ) 


( 21 ) 


because g(0) is zero (when T = TJ . In a similar fashion, it can be 
shown that 


*"(6) - - ^ (Jr 6 - jf“) + V 

P P P 




*<»T > 


( 22 ) 


P P 

After some reduction, Eq. (22) evaluated at 6 s 0 gives 


g”(0) 


(C ). • (C ) 


JtlX 


Vt 


P'o 


2o»h 4 


An examination of Eqs. (12) *(23) reveals that 


a o - 0, aj * h/T t> 


*(AC p / 2 v a h t ) ♦ (« 0 h t /T t ). ... 


(23) 


(24) 


(c) Results. 

The results obtained through the series expansion will 
now be discussed. These results will be interpreted for both first and 
second order Ehrenfest transitions. These mathematical results will 
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then ba used to give semiquantitative raault* for both clasao* of 
transitions in raprasantativa systems. 


First ordar transitions: 

A ravisitation of Eq. (9) with tha incorporation of tha 
axpansion darivad in tha previous sections gives a critical nucleus 
r c * «2T/g - -21/(8^ + s 2 0* + ...), (25) 


2r(T t /0) 


- I (AC /2v.h t ) - o ft ] 0 + ...» 


(26) 


where h^ is a negative number for an exothermic transition. 

Similarly, incorporation of Eq. (10) with tha expanded results 
gives an activation energy barrier 

AG c - (16s/3)t a /g* « (16«/3)» a /(« 1 0 + a 2 0 a + ...)*, (27) 




(16ir/3)y J (T t /0)* 

( h ^) 1 


1 - [(AC p /2v p h t ) - ep) 1 © 1 " + • 


(28) 


These equations reduce to the conventions! expressions given in 
Eq. (11) when the higher-order correction terms are neglected. Thus 
as a result of the expansion of the free energy in a Taylor series, we 
have explicit, convenient, and (in principle) experimentally accessible 
correction factors for first order homogeneous nucleation theory. 


Second-order transitions: 

In a second-order Ehrenfest transition, since the latent heat is 

zero, the conventional expressions given in Eq. (11) are physically 

meaningless because in both cases the denominator is zero. However, 

using the results derived in the preceeding section for r and AG., 

c c 



meaningful results for the critical nucleus size and critical activation 
energy for second order transitions can be determined. Thus 
evaluating Eq. (26) for the case where h^ = 0 gives a critical nucleus 


. 4! Vt 

r c (AC )0 1 * 
P 


(29) 


Similarly, evaluating Eq. (28) for the case where h^, = 0 gives 
an activation energy barrier 



(64w/3)y*u’T 

P t 

(AC )*0* * 

P 


(30) 


There is evidence that the interfacial energy Z ~ between 
phases separated by a second order transition may have a zero value at 
the equilibrium second order transition temperature. However on 

the basis of our continuing investigation we suspect that at a degree of 
undercooling 0 = T - T^, which is the point of concern in the present 
argument, J will not vanish because the curves for the total surface 
energies y(a) and 3f(&) have different slopes ay(a)/3T and ay(3)/3T 
(Fig. 3). In the case of the glass transition in polymers where 

T^ = T , a corresponds to the rubbery (liquid) state and fl to the 

glassy state. The underlying physical phenomenon is the existence of 
molecular fluctuation which give rise to the formation of a small embryo. 

Hence a mathematical formulation for applying classical nucleation 

theory to second order Ehrenfest transitions is found through the 

expansion of the Gibbs free energy in a Maclaurin series. 
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(d) Numerical results for undercooled systems. 


A semiquantitative graphical representation of the results 
is given in Figs. 4-9 using typical data for the first order liquid to 
solid transition of water and the quasi-second order glass transition of 
polystyrene. In Fig. 4 and Fig. 5 the critical parameters for a 

first order transition are plotted against the degree of supercooling 
using both the conventional and expanded approaches. These plots 
serve to give an order of magnitude estimation as to the size of the 
"correction factor" discussed earlier. Figure 6 shows a relationship 
similar to the qualitative representation of Fig. 2 in which the free 
energy is plotted against the cluster size for various degrees of 
supercooling. The peak of each curve gives the values for r and G . 

C V 

The values of the data were collected from various sources in the 

literature and are listed in Refs. 11 and 12. 

Figures 7 and 8 graphically show the size of the critical 
parameters as a function of assumed interfacial energy between the 
glassy and rubbery states of polystyrene at a fifty degree super- 
cooling. The use of an assumed interfacial energy is a consequence of 
the lack of data of this type in the literature. The values chosen, 
although somewhat arbitrary, are believed to be of the proper 

magnitude for such systems. Figure 9 shows the variation of the 

critical nucleus size with the assumed interfacial energy for various 
degrees of supercooling. The source of the data plotted for 
polystyrene is given in Refs. 13 and 14. 

As indicated above, appropriate experimental data for the 

interfacial energy = Y are not presently available. Therefore we 
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have chosen Y as an independent variable for graphical representation 
in Figs. 7, 8, and 9. Current research efforts are directed toward the 
determination of appropriate values of Y. These plots serve to provide 
guides for the design of relevant experiments to measure this interfacial 
energy. 

2. Expansion in Terms of Applied Hydrostatic Pressure. 

The formalism developed earlier and applied in the previous 
section will now be utilized to predict the effects of an applied 
hydrostatic pressure. 

(a) Predicted effects of hydrostatic pressure. 

As outlined earlier, the main point of the theory is the 

expansion of the free energy in a Maclaurin series. 

This expansion in terms of AP = P-P , where P is the initial 

o o 

pressure of the system, may be expressed as follows 

g = g(AP) = g(0) + g'(0)AP + ig"(0)(AP) a + ... (31) 

= b Q + bjAP + b 2 (AP) 2 + ... (32) 

where, for example, b 1 = g'(0) = (3g/aP)j evaluated at AP = 0. It 

should be noted that g is also a function of temperature and thus g 

evaluated at AP = 0 is not necessarily zero unless the function is also 

evaluated at the normal temperature for the transition at the initial 

pressure of the system. Alternately stated, the function is also 

evaluated at a degree of supercooling 0 equaling zero as in BC-I. 

Thus for example a system which is initially at atmospheric pressure 

(P = 1 atm.), g(0) equals zero when the function is evaluated both at 
P - P q or AP ~ 0 and at the normal transition temperature T = T^ or 0 



= 0. Since this is the case in the present discussion, b Q * g(0) in Eq. 
(32) will be zero. 

From the definition g s AG/u a and relations such as u s 

P a 

(3G /dP)--, and B = -(3lnu /3P) T for the volume and compressibility of 
a l a a l 

the a-phase, it follows that 


g'(0) 5 


( -2S) 

l 9AP'. 


(is, .-l !<!!§) 

»T v » «T 


8G 


G r -G 3 v r 


(33) 


AP=Q 


Thus at AP = 0 and 0 - 0 we have 


g'(0) = b 1 


- V P V g _ Av 

V P V P 


(34) 


where Ao = u a - u . Note that G a - G in Eq. (33) is zero for the 
pa pa 

reasons discussed previously. 

In a similar way. 


i r 


M , 0 -v = 1 b n (35) 

g (0) v & | ( 3P ) T ( 3P ^1 " ^| (v p" V a )( 3P } T ! | g( lP^ T ( 1 P t B P|| 


I I 

J L 


AP=0 


After reduction and evaluation at AP = 0 and 0 s 0 

r i 

|v ( v R ~v )P R j 

ig"C°) = b 2 = ij- (S o - *„) + -*r— *J 

L P P J 


(36) 


Summarizing the results of the expansion 

g(0) = b Q = AG| ~ 0 

v e I AP=0 


(37) 


v„-v 


g' (0) = b. 

1 V P V P 


(38) 
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(39) 


V 

b * *(- fit + ~Av) 

44 V P V & 

(b) Results for a first-order transition. 

Use can now be made of Eqs. (9) and (10) by 
incorporating the results of the free energy expansion to obtain 
expressions for the critical nucleus size and critical energy barrier for 
a first order transition occurring isothermally under the influence of an 
applied pressure. The critical nucleus size is found from Eq. (9) to be 


r 

c 


r 

c 


= 2*/g = 


-2r/(b o +b 1 AP+b 2 (AP) 2 + ...) 


- 2 * 


•4Sap 

v e 


L P P J 


Similarly the critical energy barrier is found to be 


AG 

c 


(i6ir/3)y 3 

r v B 

* ifcpAA - — Av) (AP) 
B V B 


1 2 

\ 

I 

J 


(40) 

(41) 


(42) 


Thus relatively simple expressions have been found for r and AG 

c c 

for phase transitions induced by the application of pressure. 


(c) Results for second-order transitions. 

As in the previous section, expressions for r c and A(^ 
can be found for second-order transitions which occur in the presence 
of an applied pressure. For a second-order transition 

v = v- or Av = 0 (43) 

o p 

so that the value of g reduces to 


-*AB(AP) a 


14 


(44) 



Therefore the theoretical value of r in a second-order transition, 

c 

induced by applied pressure, is found to be 


„ , 4Y (45) 

r c " A&(AP) 2 

Similarly, the theoretical value for the critical energy barrier is 




(46) 


3. Expansion in Terms of Applied Tensile Stress. 

The new approach will next be applied to the case of an 
applied tensile stress. 

(a) Mathematical development in terms of stress. 

The same type of derivation which was carried out for 
the case of applied pressure will now be developed for the application 
of a tensile stress. Although very similar in some respects to the 
previous case, there are some important differences which require a 
rather detailed presentation, rather than merely treating the stress as a 
negative pressure. Starting in the same way as for Eq. (31) we have 
g = g(o) = g(0)+g , (0)a+ig ,, (0)o 2 -t- ... (47) 

= b' + b' o + b' o 2 + ... C48) 

0 1 2 


where o = F/A o is the nominal stress and g'(o) denotes (3g/3o)p j at o 
= 0. This derivative can Ire expressed as follows 


(JLS) 


P,T 


= (IS.) 

''3F' 


P.T 


(il) 

' > 30 ; 


= A 


P.T 


( i*> 

> k 3p 


P.T 


( 49 ) 


where A is the initial cross-sectional area, 
o 
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Carrying through the required operations 


* * i 36 G a -G 3v ft 

l’W) s (!?) = Aclf) = A< i (T=t) - (iA - (-MitjS 6 


3o 3F' 

P,T 


v« I 3F 7 3F I T v« * 3F ' |'| 

P,T P | P,T P,T| j P P,T j j 0=0 

L J L J 


Use will now be made of the thermodynamic relation (3G/3F)p j - -L. 


Thus at o s 0 or F = 0 we have 


• " ( V L «> -AT 

g‘(0) = b = A — -- A — - F 1 - = — A 

1 ! ’» ] v 6 

where AL = L a -L and g(0) is zero when o = 0 or F = 0. 
p a 

In a similar fashion, it can be shown that 


-2 | . [ 3L 3L I f(L -L ) 3 V - 1 

8 CO) Cgp2)A | A l ir I ( 9r ) C ar ) p T | j T C air ) j 


= a I ( SLi 


S\ 3F P,T 


3F ^ P ,T | 


3F I 

D T* I 


. OV fi j 

'**5^ -8^ I 

P,T P,T| 


lji ( !V + iei) ( i!§) «> 

| ( 3F } l^: )( 3F^ ) p T l } 

i P,T j P P,T, o=0 


After reduction, Eq (52) evaluated at o = 0 gives 


r i r i 

I 1 1 3v ! 

_ h * _ f 1 I , L . /i J . t f AL I 

ig (0) b 2 - ( 2 |( YA ) ( YA ) | + | ( T“* )( 3F^ | 

Pi_ a Pj l P *» 


where use has been made of the relation (3L/3F) = L/YA_ where Y is 

o 


the Young’s modulus. 

Collecting the resulting coefficients. 


g(0) = b = — ! = 0 

' o Vp | o=0 


16 



(55) 


g’(0) 


t 



- ( vv 



i*"(o) = »; - ( 2 i-jc^) 

P j _ o 


f Ji v I 


L P 


) )A a 
P.Tj 


(56) 


(b) Results for a first order transition. 

Followin', the procedure outlined previously for the case 

of applied pressure, one can find the expressions for r and AG for 

c c 

the case of applied stress. 

r c = -22T/g = -2I/(b , ( +b i , 0+t/o 2 + ...) (57) 


r 

c 


■ 2 * 


MAo r 

_( v fi ) + { 2v s | ( YA ) 

P Pl « 


V-A) I 

YA> sJ [ v 8 8F p,tJ 


2*2 


}o 2 A 


(58) 


Similarly the critical energy barrier is 

AG = (16w/3)* , /g 55 ( 16if/3) y* / (b * +b 1 o+b 1 o 2 + ...)* (59) 

c o 

1 2 


AG a 
c 


ii6*/3)y»_ 


-ALAo 


i r l 


+ Cya) & J 


[AL 3 v - 1 

o 2 A 2 + !- (,/) Jo l A 2 ) 2 

P»Tj 


(60) 


v » 


Thus through the expansion of the free energy in a Maclaurin 
series one can derive expressions for the critical nucleus size and 
critical energy barrier for a system undergoing a first order phase 
transition under the influence of an applied tensile stress. 
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(c) Results for a second-order transition. 

Likewise the expressions for r and AG can be found for 

c c 

a second-order transition by adjusting the first order expressions to be 
consistent with the thermodynamic criteria of second order transitions. 

For a second order transition, 

L * L fl or AL * 0 (61) 

o p 

Thus the theoretical value of the critical nucleus size in a 
second order transition in the presence of a tensile stress is 


r 

c 


4 , v 6 

s ■ »JL i r 
o 3 (AY) 


(62) 


where AY = Y ft -Y . 

P « 

The theoretical value of the critical energy barrier is then 
found to be 


(64s/3)Jr*Y* Y* 
AG c s o- (AY) t 


(63) 


4. Investigations of Polymer Solid-Liquid Interfacial Energy. 

Due to the importance of the interfacial energy term in the 
expressions derived in the previous sections, a great deal of effort was 
put forth searching the literature for such data. As mentioned earlier, 
the interfacial energy between two "phases" in a second order transition 
is a quantity which customarily is believed to be zero at the normal 
transition point. However, as explained, the value may be non-zero at 
a degree of undercooling or under an applied pressure for example. 
Because this way of looking at this interfacial energy is somewhat novel 
and arises due to our particular interest in this situation, it is not 
surprising that mention of this quantity is not found in the literature. 
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However, what is surprising is that there is relatively little data 
available on the interfacial energy between the solid and liquid phases 
(Tg|_) in polymeric materials. As a result of this lack of data and the 
information collected in search if a value for the glass-rubber 
interfacial energy, an investigation into the problem of evaluating Yg^ 
was undertaken. 

(a) Methods of determining Yg^ for polymers. 

In this investigation, four methods were used to evaluate 
the liquid-solid interfacial energy. In method I, the interfacial energy 
was computed using the simple relationship 



( 64 ) 


where Yg is the solid surface tension and Y^ is the liquid surface 
tension. In this method Y^ was found at the melting point T m through 
the use of a large amount of available data for liquid surface tension 
and (3Y L /aT). The relationship 


Pc n 

r s * V 


( 65 ) 


where pg and p^ are the densities of the solid and liquid phases and n 
is the MacLeod exponent (generally found to be about 4), was then 
used to compute Yg. From the difference between Yg and Y^, the value 
of Yg^ is then found. The results of this method for various polymers 
are shown in Fig. 10a plotted against the molar heat of fusion. This 
type of plot is based on the work of Turnbulr ' on metals and Thomas 
and Stavely^^ on organics where similar plots for these materials yield 
straight line relationships. As a result of these linear correlations, the 
solid-liquid interfacial energy can easily be estimated for a class of 
materials with knowledge only of the molar heat of fusion, a quantity 
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which is much more readily available in the literature and obtainable 
experimentally. 

As can be seen in Fig. 10, the first method does yield a somewhat 
linear relationship. However, this method does not account for the 

f 181 

polarity of the polymer and thus is probably not very accurate. 

In method II, the polarity is considered by using a formula 
derived by Wu^®^: 


SL 




4* d r d 

r d +y d 

*S L 


p p 
# S # L 


( 66 ) 


where Yg and Y^ are the dispersive components of the solid and liquid 

P P 

surface tension and Yg and Y^ are the polar components of the solid 

(18) 

and liquid surface tension. Using values given by Wu v for the 

quantities in Eq. (66) the result for Yg^ is shown in Fig. 10a. 

The third method used to estimate Yg^ is that of Owens and 
(19) 

Wendt . Their equation is 



Y 

c 


(67) 


where Y is the critical surface tension. Again using data given by 

V 

(18) 

Wu v ’ for Y , the results are shown in Fig. 10b. It should be pointed 
c 

out that the values of Y which were used are at 20°C and not at T 

c m 

and thus some error may have been introduced due to this fact but it 
is thought that the plot provides useful information. 

The final method (IV) for estimating Yg^ is the same as that 
described by Turnbull^®^ and Thomas and Stavely^^. In this 

method, an estimate is made for the activation energy barrier and 
Eq. (11) is solved for the interfacial energy Yg^ for a certain degree of 
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( 16 ) 

undercooling 6. Using similar approximations to those of Turnbull 
and Thomas and Stavely^^, the results obtained for are shown in 
Fig. 10c. 

This fourth method appears to give the best results and is by far 
the easiest of the four methods to obtain data and estimate *SL* 
Furthermore, by refining the assumptions made in the analysis, even 
greater accuracy may be possible. 
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SECTION IV 


SUMMARY AND COMMENTS 

In this work we have developed a formalism for predicting the 
critical nucleus size and activation energy for second order Ehrenfest 
transitions. Correction terms are also derived for first order 
transitions which break from the customary approximations normally 
employed in classical nucleation theory. In completed work, the effects 
of undercoolng, hydrostatic pressure, and externally applied tensile 
stress have been rigorously investigated. In the case of temperature 
undercooling, numerical values for the critical quantities were obtained 
using a somewhat arbitrary range of values for 1 the interfacial energy. 
No such attempt was made for the cases of applied pressure and stress 
due to the uncertainty of valu*f for this term. Later work on this 
problem may have given a possible solution and method of estimating the 
interfacial energy between two "phases" undergoing a second order 
transition. This method concerns the difference in slopes of the Jf vs. 
T curves for the different phases as described in this report. However 
due to the lack of t vs. P and I vs. o data, such an estimate is still 
not possible for the other two situations described herein. 

Work is now being done in which a Claussius-Clapeyron relation is 
used to attempt to find I vs. P data and thereby reasonable values for 
the interfacial energy in the case of applied hydrostatic pressure. 

Work is also continuing in which the effects of an applied electric 
and magnetic field are considered. It is hoped that such a situation 
will be more conducive to experimental investigation and thereby yield 


24 



experimental evidence for the ascertions made both in completed work 
and work still to be undertaken. 

A number of papers and presentations have resulted from work 
described in this report, and our efforts are continuing in anticipation 
of many more fruitful developments as the research proceeds. 
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FIG. 1. For a phase transition a ■» gto occur at a degree of 

supercooling 0, G ft must be less than G where G is the molar Gibbs' free energy. 
P a 

A# = 4/3 7rr3 g + 4 irr^y 



A $ c = energy activation barrier = 167rx 3 /3g 2 
r c = critical nucleus radius = -2//g 


FIG. 2. The two opposing energy terms which give rise to a critical 

radius r and a critical energy barrier AG for nucleation. 
c c 



TEMPERATURE T 


FIG. 3. Schematic representation of a method to estimate the effective 

interfacial energy for an undercooled system if Z = Z . 

y * 


CRITICAL RADIUS r c (ANGSTROMS) 
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FIG. 4. The dependence of critical nucleus size on the degree of 
supercooling for the transition of liquid water to ice for both the 
conventional and expanded approaches. 
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FIG. 5. The dependence of the free energy barrier on the degree of 
supercooling for the transition of liquid water to ice for both the 
conventional and expanded approaches. 
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FIG. 6. The free energy as a function of cluster size for various 
degrees of supercooling for the transition of liquid water to ice for both 
the conventional (solid curve) and expanded (dashed curve) 
approaches. 
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ASSUMED INTERFACIAL ENERGY /(10‘ 5 J/m* ) 

FIG. 7. The critical nucleus size as a function of assumed interfacial 
energy for the glass transition of polystyrene at a supercooling of 50 
deg. 
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FIG. 8. The free energy barrier at a function of assumed interfacial 
energy for the glass transition of polystyrene at a supercooling of 50 
deg. 
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FIG. 9. Tht critical nucleus size as a function of assumed interfacial 
energy of various degrees of supercooling for the glass transition of 
polystyrene. 
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Predictions of nuclcation theory applied to Ehrenfest thermo- 
dynamic transitions : II. The effects of pressure and stress. 

K. W. Campbell and R. E. Barker, Jr. 

Department of Materials Science, University of Virginia, Char- 
lottesville, Virginia 22901 
(Received 

Abstract 

I 

This paper is a sequel to an earlier one, BC-I (J. Appl. Phys . 56, 2386 

(1984)), on the applicability of classical nucleation theory to second- 

order transitions in the Ehrenfest sense. In each case the approach was 

to obtain the critical size r and energy barrier AG for the growth of 

c c 

a nucleus of B~pha.se in ati a-phase matrix by a Maclaurin series expan- 
sion of the f ree-energy-density g=(G fl -G )/u Q as a function of 0 (in 

P o p 

BC-I) and of AP and Ao in this paper where 0=(T-T t ) is the degree of 
undercooling (and AP and Ao are analogous terms for the hydrostatic pres- 
sure shift and tensile stress shift away from the equilibrium transi- 
tion. The expansion coefficients were determined by the use of thermo- 
dynamic relationships. For second- order transitions, r^=4I\)pT t /AC^0 a , 
r c =4I/AB(Ap) J , and r^AYY^Y^/AYCAo)* , respectively for the three cases. 
The terms AC^, AB, and AY denote the differences in heat capacity, 
compressibility, and Young's modulus, e.g., AY=Yp-Y o . The interfacial 

energy H a is denoted by T. The activation energy barriers for the 
ap 

cases developed in this paper were A G £ =( 16ir/3) */ (A3) 1 (Ap) * and 

AG =(64u/3)Y 3 Y a Y fl a /(AY) a (Ao) fc . More complicated expressions are given 

c a p 

in the paper for the r £ and AG c for first-order transitions. In the 
long run these expressions may prove more useful than the ones for 
second-order because of the modifications expressions for the kinetics 
of transformations. 


1 


where V § it the volume of the p nucleus, v 0 is the molsr 
volume and 4* is the interfsciil tree of p. If the growing 
nucleus has more than one type of interface, then the term 

y* would imply a summation over the relevant areas. 
In the present paper, an extra term, corresponding to the 
elastic strain energy when v a is not equal to v 0 in solid 
phases, will be neglected. 

Thus, according to nucleation theory, local density 
fluctuations will lead to the formation of a small nucleus of/9 
{dune within an a-phase matrix if the temperature is less 
than the normal transition temperature T, for or where 

P is the stable phase when T<T, (Fig. 1). The total interfa* 
dal energy y will oppose the enlargement of the nu- 
cleus, and s bulk-free energy gV 0 will encourage growth 
(Fig. 2). The term 



FIG. 2. The two opposes energy terms which give rise to a critical radius r € 
and a critical energy barrier A 9 C for nucleation. 


g-(G 0 -G a )/v 0 (3) 

represents the free energy change per unit volume of trans- 
formed material. Equation (2) now becomes 

A&=gV 0 + yA 0 . (4) 

C. The customary approximations 

In the standard treatments, AG = AH — TAS is ap- 
proximated, at a degree of undercooling 

8 =T,-T (5) 

by AH = AH, &adAS = AH,/T, which may be seen to give 
AG = AS,9, with S, - AH { /T t (6) 

or 

g = r, ft with s, = AS,/v 0 . (7) 

The standard treatments next proceed by taking d{AG)/ 
dr = 0, where r is the radius of the spherical nucleus for 
which 

V 0 = (4/3)irr 3 and A 0 = 4flT 2 . (8) 

This gives the size of the critical nucleus r c to be 

r c = - 2y/g, (9) 

For the stable a~*fi transition, g will be negative (g = 0 for 
the equilibrium transition). Substitution of Eq. (9) back into 
Eq. (4) gives the usual type of expression for the “activation 



FIG. 1. Fort phase munition to occur at a degree of supercooling#, 
G f must be less than G„ where G i* the molar Gibb*' free energy. 


barrier” that must be overcome before the continued growth 
of the P nucleus can take place. 

A9 e = (1 (nrtltf/f = 16.76 f/g 2 . (10) 

This barrier, and the critical size r e , are of course both de- 
pendent on the degree of undercooling, through Eq. (7), 
therefore 


r e = 


— 2y/s,8 — 


-2 yjT,/9) 

h, 


( 11 ) 


and 


A& c = {16tr/3)y 3 /st0 2 = 


\\<ytrn\r 3 {T,/d ) 2 

( A ,) 2 


( 12 ) 


where h, = AH, /v 0 is the latent heat per unit volume for the 
transformation at standard pressure. 


NEW THEORETICAL APPROACH 

A description will now be given for a mathematical for- 
mulation of nucleation theory in which the customary ap- 
proximations discussed earlier are avoided. The result is a 
number of expressions which are unobtainable using the 
conventional methods. 

A. Proposal of a more rigorous a ppr oa ch 

The twelve equations above summarize an elementary 
view of homogeneous nucleation theory for the special case 
of spherical particles. The need for a more rigorous formal- 
ism in which the customary approximations are not applied 
is evident. In this new theory, such a formalism is developed 
through the assumption that A if and therefore g can be 
expanded in a Maclaurin’s series in 8 , the degree of under- 
cooling. This series is then reduced and manipulated 
through the use of Maxwell’s relations and other thermody- 
namic operations. As a result of this rigorous treatment, the 
following features are developed: (1) nonlinear correction 
terms result for first-order thermodynamic transitions, (2) 
the critical nucleus size and barrier height for second-order 
transitions (in the Ehrenfest sense), and (3) the possibility of a 
formal extension erf 1 nucleation theory to even higher-order 
transitions. Since the new theory brings in the higher-order 
effects as analytic thermodynamic terms (rather than merely 
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mathematical or empirical approximations) it is a more flexi- 
ble theory for predicting other effects, such as the influence 
of changes in pressure, stress, electric field, etc. The new 
theory, for the same reason, should also be more amenable to 
tk» of ideas from uui to 

compositional and kinetic extensions such as Couchman’s* 
T„ theory and the time- temperature-transformation con- 
cepts of Enns and Gfllham. 9 


B. Central teaturee of the new approach 

As outlined in the previous section it is possible to con- 
struct a useful, ami apparently overlooked, formalism by 
what in its barest form is the expedient of expanding the 
effective free energy of the transformed nucleus in a Ma- 
claurin series in 6 (the degree of undercooling), and then 
utilizing thermodynamical relations to evaluate the expan- 
sion coefficients. For example, if we take 

g-glfi) -«<0) +fW + (l/2)g'(0)d 2 + ... 

— a 0 + a x 6 + a^ 2 + ... (13) 

then g'(O), which means {dg/dO ) p at 6 = 0, can be found as 
follows. 



(14) 

but 


G fi -G a =H 0 -H a -T[S 0 -S a ), 

(15) 

so that 


g&{Gp — G a )/vp = h — Ts, 

(16) 

where 


h ={H 0 - H a )/v 0 and s={S p - S a )/v p . 

(17) 

Carrying through the operations required, 





(18) 

— — ■ {Sff ~ S a) ~ 8 a p< 

(19) 


where use has been made of the thermodynamic relations 
S— — {dG/8T) p and thermal expansivity = a = \dv/ 
dT),/v. 


C. Calculationc of the coeffi ci ents 

Further reduction of Eq. ( 18) gives 



and thus, using Eq. (13) at T— T t or Q — 0, we obtain 






( 21 ) 


because g(0) is zero (when T = T t ). In a similar fashion, it can 
be shown that 


After some reduction, Eq. (22) evaluated at $ ** 0 gives 


8'{0) 


(C,) 0 -(C,) a . 2 a,h t 


Vf T, T t 

An examination of Eqs. (12)— (23) reveals that 
0 q = o, a, = A/r,, 

O 2 — — {AC p /2v fi h t ) + (ctph,/T,),... 


(23) 


(24) 


RESULTS 

The results obtained through the series expansion will 
now be discussed. These results will be interpreted for both 
first- and second-order Ehrenfest transitions. These math- 
ematical results will then be used to give semiquantitadve 
results for both classes of transitions in representative sys- 
tems. 


A First-order transitions 

A revisitation of Eq. (9) with the incorporation of the 
expansion derived in the previous sections gives a critical 
nucleus 

r c = - 2 y/g = - 2y/{a t 6 + a 2 0 2 + ...), (25) 


r e = - 


2 riT,/e) 


1 


(26) 


h, l-[{AC p /2v / ,h l )-a /l ]0+... 

where h , is a negative number for an exothermic transition. 

Similarly, incorporation of Eq. (10) with the expanded 
results gives an activation energy barrier 

A& e « (16ir/3)r Vg 2 = (16jr/3)y 3 /{a l 0 + a 2 0 2 + ...) 2 ,(27) 
( 16ff/3)y 3 (T,/0) 2 




( A ,) 2 


X . (28) 

1 -[{AC p /2vph,)-a fi ] 2 0 2 + ... ' ' 

These equations reduce to the conventional expressions 
given in Eqs. (11) and (12) when the higher-order correction 
terms are neglected. Thus as a result of the expansion of the 
free energy in a Taylor series, we have explicit, convenient, 
and (in principle) experimentally accessible correction fac- 
tors for first order homogeneous nucleation theory. 


B. Socond-ordor transition* 

In a second-order Ehrenfest transition, since the latent 
heat is zero, the conventional expressions given in Eqs. (11) 
and ( 12) are physically meaningless because in both cases the 
denominator is zero. However, using the results derived in 
the proceeding section for r c and A9 e , meaningful results 
for the critical nucleus size and critical activation energy for 
second-order transitions can be determined. Thus evaluat- 
ing Eq. (26) for the case where h,~ 0 gives a critical nucleus 
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_ _ +rv* T < 

Similarly, evaluating Eq. (28) for the earn where k, 
gives an activation energy barrier 


(29) 
= 0 


_ PWjg? 

• (JC,)’»* 11 

There is evidence that the interfacial energy y — y^ 


between phases separated by a second-order transition may 
have a zero value at the equilibrium second-order transition 
temperature. 10 However mi the basis of our continuing in- 
vestigation we suspect that at a degree of undercooling 
0 = T— T„ which is the point of concern in the present 
argument, y will not vanish because the curves for the total 
surface energies y{a ) and yi p ) have different slopes dy\a)/dT 
and dy{P)/dT. In the case of the glass transition in polymers 
where T, = T t , a corresponds to the rubbery (liquid) state 
and P to the glassy state. The underlying physical pheno- 
menon is the existence of molecular fluctuation which give 
rise to the formation of a small embryo. 

Hence a mathematical formulation for applying classi- 
cal nucleation theory to second-order Ehrenfest transitions 
is found through the expansion of the Gibbs’ free energy in a 
Maclaurin series. 


w* NWnfnc« mum 

A semiquantitative graphical representation of the re- 
sults is given in Figs. 3-8 using typical data for the first-order 
liquid to solid transition of water and the quasi-second-order 
glass transition of polysytrene. In Fig. 3 and Fig. 4 the criti- 
cal parameters for a first-order transition are plotted against 
the degree of supercooling using both the conventional and 
expanded approaches. These plots serve to give an order of 
magnitude estimation as to the size of the “correction fac- 
tor” discussed earlier. Figure 5 shows a relationship similar 



SUPERCOOLING 8 (DEG) 



SUPERCOOLING 8 (DEG) 

FIG. 4. The dependence of the free energy barrier on the degree of super- 
cooling for the tranrition of lk|uid water to ice for both the conventional and 
expanded approaches. 


to the qualitative representation of Fig. 2 in which the free 
energy is plotted against the cluster size for various degrees 
of supercooling. The peak of each curve gives the values for 
r c and G e . The values of the data were collected from various 
sources in the literature and are listed in Refs. 1 1 and 12. 

Figures 6 and 7 graphically show the size of the critical 
parameters as a function of assumed interfacial energy 
between the glassy and rubbery states of polystyrene at a fifty 
degree supercooling. The use of an assumed interfacial ener- 
gy is a consequence of the lack of data of this type in the 
literature. The values chosen, although somewhat arbitrary, 
are believed to be of the proper magnitude for such systems. 
Figure 8 shows the variation of the critic*! nucleus size with 



FIO. 3. The de p e n de nc e of critical nucleus tiie on the degree of supercool- FIG. 5. The free energy as a function of cluster iize for various degree* of 

lag for the transition of liquid water to ice for both the conventional and supercooling for the transition of liquid water to ice for both the convention- 

expanded approaches. al (solid curve) and expanded (dashed curve) approaches. 
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ASSUMED INTERFACIAL ENERGY y(10‘* J/m* ) 


FIG. 6. The critical nucleu* size m a function of aaaumed interfacial energy 
for the gleet transition of polystyrene at a supercooling of 30 deg. 

the assumed interfacial energy for various degrees of super* 
coding. The source of the data plotted for polystyrene is 
given in Refs. 13 and 14. 

As indicated above, ap pro p riate experimental data for 
the interfacial energy y^ == y are not presently available. 
Therefore we have chosen y as an independent variable for 
graphical representation in Figs. 6, 7, and 8. Current re* 
search efforts are directed toward the determination of ap- 
propriate values of y. These plots serve to provide guides for 
the design of relevant experiments to measure this interfacial 
energy. 


CONCLUSION 

In this paper we have developed a formalism for pre- 
dicting the critical nucleus size and the activation energy 
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ASSUMED INTERFACIAL ENERGY y(10'° J/m*) 


FIG. 7. The free energy barrier at a foncuon of aisumed interfacial energy 
for the glee* transition of polystyrene at a supercooling of 30 deg. 



ASSUMED INTERFACIAL ENERGY y<10** J/m*) 

FIG. 8. The critical nucleus si» at * function of assumed interfedal energy 
of various degree* of supercooling for the fleet transition of polystyrene. 


barrier for second -order Ehrenfest transitions. Correction 
terms are also derived for first-order transitions which break 
from the customary approximations normally employed in 
classical nucleation theory. A semiquantitative analysis of 
this new approach was that presented to give physical 
sense to the theoretical relations. Thus as a result of expand- 
ing the Gibbs’ free energy in a Maclaurin series, a mathemat- 
ical model is developed for applying classical nucleation the- 
ory to second-order Ehrenfest transitions. 
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Abstract 


f 


This paper is a sequel to an earlier one, BC-I (J. Appl. Phys . 56, 2386 
(1984)), on the applicability of classical nucleation theory to second - 
order transitions in the Ehrenfest sense. In each case the approach was 
to obtain the critical size r £ and energy barrier AG^ for the growth of 
a nucleus of &-phase in an ct-phase matrix by a Maclaurin series expan- 
sion of the free-energy-density g=(G 0 -G )/u Q as a function of 0 (in 

pa p 

BC-I) and of AP and Ao in this paper where 0=(T-T ) is the degree of 
undercooling [and AP and Ao are analogous terms for the hydrostatic pres- 
sure shift and tensile stress shift away from the equilibrium transi- 
tion. The expansion coefficients were determined by the use of thermo- 
dynamic relationships. For second* order transitions, i v LC ^Q 2 , 

r c =42f/AE(Ap) z , and r c =4*Y a Yp/AY(Ao ) 2 , respectively for the three cases. 
The terms AC , AP, and AY denote the differences in heat capacity, 

p 

compressibility, and Young's modulus, e.g., AY=Y^-Y a . The interfacial 

energy Z _ is denoted by Z. The activation energy barriers for the 

op 

cases developed in this paper were A G^=( 16 tt/ 3) Jf J / (Af3 ) 7 (Ap) * and 

AG =(64it/3)r J Y 1 Y e l /(AY) 2 (Ao)* . More complicated expressions are given 
c a P 

in the paper for the r^_ and AG c for first* order transitions. In the 
long run these expressions may prove more useful than the ones for 
second-order because of the modifications expressions for the kinetics 
of transformations. 
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Introduce ion 


In a previous paper (BC-I) 1 on this general subject, a formalism 
was developed for the application of classical nucleation theory to 
firsthand second-order thermodynamic transitions in the Ehrenfest 
sense. This was accomplished through the expansion of the Gibbs' free 
energy in a Maclaurin series in the degree of supercooling 0. 3y using 
this approach, significantly improved expressions were obtained for the 
size r^ and the energy barrier A of the critical nucleus for first- 
order transitions. Moreover, physically meaningful expressions also 
were obtained for r c and A for a system undergoing a second- order 
phase transition with the assumption that, in such a transition (e.g., 
order-disorder , normal-superconducting, and (approximately) the glass 
transition in polymers), the new phase does not occur spontaneously 
throughout the matrix but rather begins locally as nuclei which grow 
until the transformation cr*3 is complete. For this approach, the 
existence of an effective interfacial energy 7=3^ between the two 
"phases" (or states of aggregation) is essential. However in the usual 
theories of second- order transitions X is believed to be zero at the 
equilibrium transition temperature T . 1 For the purposes of the present 
analysis we postulate that the type of transformation of interest (a+fl) 
occurs under non-equilibrium (supercooled) conditions with a small but 
nonvanishing X. Furthermore, by the application of pressure, tensile 
stress, or electric, or magnetic fields of sufficient magnitude, the 
free energy curves and G^ (for certain first-order transitions) which 
cross at more that one point may possibly be shifted to produce a common 
tangent, thus simulating a second- order transition. These situations 
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(supercooling and shifted free energy) thus lead to the possibility of 
nucleation and growth processes for quasi-second-order transitions. 

A formalism similar to that utilized in BC-1, wherein the excess 

free energy density g=(G c -G )/u c was expressed in terms of the degree of 

pop 

undercooling, will now be developed in terms of applied pressure p and 
tensile stress o. The results for the application of electric and 
magnetic fields will be given in a subsequent presentation. The present 
derivation, in terms of p and o, will be preceded by a brief outline of 
certain relevant features of classical nucleation theory (ChT), with 
references to appropriate literature sources for more detailed 
descriptions of the ideas and facts, especially about nucleation 
phenomena involved with polymeric materials. 

The paper concludes with a discussion of the results that are 
obtained and the implications of these results when applied to real 
systems . 

Review of Classical Nucleation Theor y 

Before discussing the detailed derivation of the expressions of 
interest in this work, a short review of CNT as it applies to this case 
is in order. For a more complete description of classical nucleation 
theory and phase transitions especially as applied to polymeric 
materials see references 1, 3-9. 

A. The thermodynamic basis of nucleation theory 

As is widely known, when a pure system of molecules which exists as 
a stable phase a at a given set of intensive variables (Tj»pj. etc.), is 
subjected to new intensive conditions (T *p * etc.), where c is no 
longer the only stable phase but where one or more new phases are 
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stable, then it is possible thnt the expected transformation a**B will 
occur only slowly, if at all, within a finite observation period. 
According to the precepts of nucleation theory, the. reason for this 
inhibition of the expected transformation is that molecular fluctuations 
and diffusive motion must occur which create a small particle of B 
within the o milieu (but see Ref. 9) and that there is an interface 
between the a and B regions in which the molecules have an excess free 
energy 1=1^ per unit area. Thus the true Gibbs free energy change 
corresponding to the transformation of a small amount of a into a small 
nucleus of B is not merely the so-called bulk thermodynamic value AG 
(energy per mole); it is instead 

AG=v b C4G/ U|1 ) + A (i 6 I i6 Cl) 

where V„ is the volume of the B nucleus, o a is the molar volume and A a 
P pop 

is the interfacial area of B- If the growing nucleus has more than one 
type of interface, then the term would imply a summation over the 
relevant areas. In the present paper, an extra term, corresponding to 
the elastic strain energy when is not equal to in solid phases, 
will be neglected. 

Thus, according to nucleation theory, a small nucleus of B phase 
will form within an o-matrix if the free energy curve of the B phase is 
lower than the free energy curve of the a-phase. This can be 
accomplished in several ways, for example by supercooling below the 
normal transition temperature T (Fig. la), or by the application of an 
external pressure (Fig. lb), or the application of a tensile stress 
(Fig. lc). 


A 



will 


Under such conditions the total interfacial energy A C T e 

ctp np 

oppose the enlargement of the nucleus and a bulk free energy gV^ will 
encourage growth (Fig. 2). The term 

* 5 <V C . ,/U e (2) 

represents the free energy change per unit volume of the transformed 
material. Equation (1) now becomes 

lG= gVp+IAg . (3) 

B. Results for s spherical nucleus 

In much but not all of the literature the nuclei of p-phase are 
assumed to be spherical for simplicity. This assumption will be used in 
this paper since we are presently more concerned with the qualitative 
features of the phenomenon, and its very existence, than with detailed 
quantitative refinements, related to nucleus shape, which we will 
address in subsequent presentations. 

In order to determine the critical parameters and A shown in 
Fig. 2., the free energy ia differentiated with respect to the radius r 
and the result is set equal to z.ero. 

6(AG)/6r=0 at r=r c * (4) 

The resulting expression for r is termed the critical radius, i.e., 
the radius that an embryo must exceed in order for its growth to be 
favorable. For a spherical nucleus where 

Vp“(^/3)TTr J and Ag-4Trr* (5) 
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the size of the critical nucleus is found to be 

r «-2»/g (6) 

c 

This value, substituted into Eq. (A), gives the critical activation 
energy barrier, 

AG =(16Tr/3)y*/* J , (7) 

c 

that is the energy barrier which inhibits the continued growth of the 9 
phase . 

New Ext ensi o n of the Theo ry 

The preceeding basic precepts of nucleation theory can be applied 
* 

to both first* and second* order thermodynamic transitions under the 
influence of applied hydrostatic pressure and more complicated stress 
f ields . 

A. Basic concepts of the new approach 

In this new approach a formalism for predicting the effects of 

applied pressure and applied stress on r and AG for a phase transition 

c c 

of any order is developed through the assumption that AG and therefore g 
can be expanded in a Maclaurin series. The coefficients in the series 

are then obtained through the use of Maxwell's relations and other 

thermodynamic operations. Relatively simple expressions result, in- 
volving experimentally accessible parameters, which describe the effect 
of applied pressure and stress on the critical quantities AG^. and r . 
More importantly these expressions yield physically meaningful 

(non-zero) results for A arid r & when the derived results are reduced 
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in accord with the theory of second order transitions. These two main 
aspects of the theory will now be described in detail, first for the 
case of applied pressure and then for the case of applied stress. 

B. Predicted effects of hydrostatic pressure 

As outlined earlier, the main point of the theory is the expansion 
of the free energy in a Maclaurin series. 

This expansion in terms of AP=P-P o> where P q is the initial pres- 
sure of the system, may be expressed as follows 

g = g(Ap) = g(0)+g'(0)Ap+J-g"(0)(AP) 2 +... (e) 

= b 0 +h,Ap +h 2 (Ap) J +... ( 9 ) 

where, for example, bj=g' (0) = (3S/3p)^, evaluated at AP=Q. It should be 

noted that g is also a function of temperature and thus g evaluated at 

AP=0 is not necessarily zero unless the function ic also evaluated at 

the normal temperature for the transition at the Initial pressure of the 

system. Alternately stated, the function is also evaluated at a degree 

of supercooling 0 equaling zero as in BC-I. 1 Thus for example a system 

which is initially at atmospheric pressure (P o =Tatm.), g(0) equals zero 

when the function is evaluated both at P=P or AP^O and at the normal 

o 

transition temperature T=T° or 0=0. Since chis is the case in the 
present discussion, b o =g(0) in Eq. (9) will be zero. 

From the definition g=AC/0 D and relations such as u =(3G /3P) , and 

P Bui 

g =-(31ri0 /3P)_ for the volume and compressibility of the a-phase, it 
a a T 

follows that 
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( 10 ) 


Ap=0 


Thus at AP=0 and 0=0 we have 


g'(o)=b, = 


V / 3 * V a Av 


( 11 ) 


where Au=u_-0 . Note that G/vG ln E< 1 ( 10 ) is zero for the 

p ° f*j a 

reasons discussed previously. 

In a similar way. 


g"(0 )= - 



1 




fev, 


pi 


( *T»\ap 



( 12 ) 


ap=o 


After reduction and evaluation at AP=0 and 0-0 


5-g"(0)= b 2 = 


1 

2 


v (v 0 ~v )B n 

— </3 -&>) + ■ ^ 0 /3 


( 13 ) 


P 


P 
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Summarising the results of the expansion 


% 


g(0) = b o = Ag 

v /3 


=0 


AP=0 


( 14 ) 
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c. 


Results for a first-order transition 


♦ 

Use can now be made of Eqs.(6) and (7) by incorporating the results 
of the free energy expansion to obtain expressions for the critical 
nucleus size and critical energy barrier for a first-order transition 
occurring isothermally under the influence of an applied pressure. The 
critical nucleus size is found from Eq.(6) to be 


«c= 2 y /g = ~2y/(b 0 +b 1 AP + b 2 (AP) 2 +... ) (17) 






(AP f + 


( 18 ) 


Similarly the critical energy barrier is found to be 


A-5[ r = 


( 1677 -/ 3 )/ 


Av 


Ap- V^Ae--£Av) (AP) 

2\ VyS P V/3 / 


2 Ci9) 
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Thus relatively simple express Loris have been found for and AG, 
for phase transitions induced by the application of pressure. 

D. Results for second-order transitions 

As in the previous section, expressions for r^ and A can be found 
for second-order transitions which occur in the presence of an applied 
pressure. For a second-order transition 

V -\/n OR Av=0 (20) 

a p 

so that the value of g reduces to 

_1A/3(Ap) 2 (21) 


Therefore the theoretical value of r in a second-order transition, 

c 

induced by applied pressure, is found to be 


A 


r c 


7 


/V 


A/3(Ap) 


( 22 ) 


Similarly, the theoretical value for the critical energy barrier is 


ah 


/v 


c 


( 6477/ 3) y 3 

(A,8) 2 (Apf 


( 23 ) 


E. Mathematical development in terms of stress 

The same type of derivation which was carried out for the case of 
applied pressure will now* be developed for the application of a tensile 
stress. Although very similar in some respects to the previous case, 
there are some important differences which require a rather detailed 
presentation, rather than merely treating the stress as a negative 
pressure. Starting in the same way as for Eq. (8) we have 


g = g( cr) = g( 0 ) + g '(O) cr + 1 g "( 0 ) cr 2 + .. . (24) 

. -bo + t) l cr + b’ 2 <x 2 + ... (25) 


where o=F/A q is the nominal stress and g' (o) denotes (3g/3o)p ^ at a-0. 
This derivative can be expressed as follows 


(§?) = eff) = A (ff) 

3 P,T 3F P,T a P,T F P.T 


( 26 ) 
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where A is the initial cross-sectional area, 
o 

Carrying through the required operations 


nm-* 


j p 


d F /pj V F /p,tj 


' Go_G a \/ dv a 


v /3 


,6F 


*.T 


(27) : 


cr=0 


Use will now be made of the thermodynamic relation (3G/3F) =-L. 

* 9 * 

Thus at o=0 or F=0 we have 


g'(0)= h', - A 


(L /T L a ' 
V /3 


-Al 


( 28 ) 


where AL=I. fi -L and g(0) is zero when o=0 or F=0. 
p a 

In a similar fashion, it can be shown that 


c"(0) = 



=a 2 ( 


cr=0 



(L o~ L /3 ) /c^A 

v* IdF / 


/P.TJ 


f bg\ 

V /3 W, 


P.T 


6F 


P.T 


73 




(29) 


cr=0 
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After reduction, Eq (29) evaluated at o-O gives 



where use has been made of the relation (3L/3F)=L/YA where Y is the 

o 

Young's modulus. 

Collecting the resulting coefficients, 


g{°)-b'o=f 


AG 


cr-0 


= 0 


g'(0)=b' = A 




( 31 ) 


-Al 


A 


( 32 ) 



F. Results for a first order transition 

Following the procedure outlined in section C, for the case of 

applied pressure, one can find the expressions for r and AG for the 

c c 

case of applied stress. 
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r c = - 2 7/g = - 2 //(b 0 + b, cr + bgcr 2 +. . .) 


(34) 


lr 


Al_Acr\ 


) 


2v, 




( 35 ) 


2 a2 


o’ A 


Similarly the critical energy barrier is 


A# c = (l67r/3) 7 3 /g = (l67r/3) 7 3 /(b 0 + b 1 cr +b 2 cr 2 -!-...) 2 (35 


) 




(167T/3)7 


(37 


( - ALAcr 1 

1 

^jj 

2 a 2 

cr A 4. 

Al_ 


cr 2 A 2 

v /3 ’ 2v £ 

\) a \h 


i 

V dF /P»T 

, 


Thus through the expansion of the free energy in a Maclaurin series 
one can derive expressions for the critical nucleus size and critical 
energy barrier for a system undergoing a first-order phase transition 
under the influence of an applied tensile stress 
G. Results for a second-order transition 


Likewise the expressions for r and AG c can be found for a 


second-order transition by adjusting the first- order expressions to be 
consistent with the thermodynamic criteria of second-order transitions. 


IS 


For a second-order transition, 


L a = Lo °R AL = 0 <m) 


Thus the theoretical value of the critical nucleus size in a 
second-order transition in the presence of a tensile stress is 


iv » 


4 r YqY/3 


09 } 


cr z ( AY) 


where AY=Y a -Y . 

P . ° 

The theoretical value of the critical energy barrier is then found 


to be 

Av*ic 


(64 77-/3) y 3 Y|Y| 
<r 4 (AY) z 


(* 0 ) 


Discussion and Conclusions 

The mathematical basis for the application of classical nucleation 

theory to second-order thermodynamic transitions has been demonstrated. 

However, a number of points should be addressed in connection with the 

derived results. The expressions for r and AG in both the case of 

c c 

applied pressure and stress contain an interfacial energy terra I. In 
the case of an equi 1 ibrium second-order transition this term vanishes.* 
However, our contention is that under non-equilibrium conditions, such 
as described here and in Ref. 1, there is a small interfacial energy 
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difference associated with the region between the a and ^-phases, for 
example between aggregations of different degrees of order. The 
question then arises whether such a transition should still be 
considered to be second-order in nature. It is our belief that the 
order classification of the transition should remain the same due to the 
fact that the transition itself is unchanged. However the kinetics are 
likely to be considerably altered. 

Our analysis has shown that thermodynamic concepts do not appear to 
preclude the nucleation and growth of a new phase in a second- order 
transition in Ehrenfest's classification. In fact reasonable 
expressions for the size and energy barrier in terms of measurable 
quantities result. However it may be that the likelihood is small for 
actually finding a set of conditions for which the predicted nucleation 
and growth can be observed. Even if this turns out to be the case, the 
results that we have derived suggest quite significant effects on the 
kinetics of first-order transitions. 

It should also be noted that with a slight modification the 
equations can be applied to the pseudo-second-order glass transition in 
polymers. Compared with the so-called genuine second-order transitions, 
the nature of the glass transition is anomalous with regard to the 
direction of tne discontinuous jumps of various thermodynamic 
quantities, such as the heat capacity and isothermal compressibility, 
which appear in the equations for r £ and A G c - The changes in Cp and & 
at T have signs opposite to the changes in these and similar quantites 
at the transition temperature for genuine second-order transition. 1 " 

In this paper a formalism is developed for the application of 
nucleation theory to first- and second-order phase transitions in the 
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Ehrenfest sense. As a result of this approach expressions for the 
critical nucleus size and critical energy barrier were obtained for the 
influence of applied pressure and stress on the transitions. There are 
also implications of important effects on the kinetics of first-order 
transitions . 
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F i gure Captions 


Fig. 1: 


For a phase transition o ), (J to occur G c must be less than G 

P o 

where G is the molar Gibbs’ free energy. This can be accom- 
plished by (a) supercooling below the normal transition 
temperature T , (b) the application of an external pressure 

AP, or (O the application of a tensile stress c 


Fig. 2: The two opposing energy terms which give rise to a critical 

radius r and a critical energy barrier AG for nucleation and 
c c 

growth . 
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